Molecular dynamics-based identification of binding pathways and two distinct high-affinity sites for succinate in succinate receptor 1/GPR91.

SUCNR1 is an auto- and paracrine sensor of the metabolic stress signal succinate. Using unsupervised molecular dynamics (MD) simulations (170.400 ns) and mutagenesis across human, mouse, and rat SUCNR1, we characterize how a five-arginine motif around the extracellular pole of TM-VI determines the initial capture of succinate in the extracellular vestibule (ECV) to either stay or move down to the orthosteric site. Metadynamics demonstrate low-energy succinate binding in both sites, with an energy barrier corresponding to an intermediate stage during which succinate, with an associated water cluster, unlocks the hydrogen-bond-stabilized conformationally constrained extracellular loop (ECL)-2b. Importantly, simultaneous binding of two succinate molecules through either a "sequential" or "bypassing" mode is a frequent endpoint. The mono-carboxylate NF-56-EJ40 antagonist enters SUCNR1 between TM-I and -II and does not unlock ECL-2b. It is proposed that occupancy of both high-affinity sites is required for selective activation of SUCNR1 by high local succinate concentrations.

Comparison of force fields to study the zinc-finger containing protein NPL4, a target for disulfiram in cancer therapy.

Molecular dynamics (MD) simulations are a powerful approach to studying the structure and dynamics of proteins related to health and disease. Advances in the MD field allow modeling proteins with high accuracy. However, modeling metal ions and their interactions with proteins is still challenging. NPL4 is a zinc-binding protein and works as a cofactor for p97 to regulate protein homeostasis. NPL4 is of biomedical importance and has been proposed as the target of disulfiram, a drug recently repurposed for cancer treatment. Experimental studies proposed that the disulfiram metabolites, bis-(diethyldithiocarbamate)­copper and cupric ions, induce NPL4 misfolding and aggregation. However, the molecular details of their interactions with NPL4 and consequent structural effects are still elusive. Here, biomolecular simulations can help to shed light on the related structural details. To apply MD simulations to NPL4 and its interaction with copper the first important step is identifying a suitable force field to describe the protein in its zinc-bound states. We examined different sets of non-bonded parameters because we want to study the misfolding mechanism and cannot rule out that the zinc may detach from the protein during the process and copper replaces it. We investigated the force-field ability to model the coordination geometry of the metal ions by comparing the results from MD simulations with optimized geometries from quantum mechanics (QM) calculations using model systems of NPL4. Furthermore, we investigated the performance of a force field including bonded parameters to treat copper ions in NPL4 that we obtained based on QM calculations.

TRAP1 S-nitrosylation as a model of population-shift mechanism to study the effects of nitric oxide on redox-sensitive oncoproteins.

S-nitrosylation is a post-translational modification in which nitric oxide (NO) binds to the thiol group of cysteine, generating an S-nitrosothiol (SNO) adduct. S-nitrosylation has different physiological roles, and its alteration has also been linked to a growing list of pathologies, including cancer. SNO can affect the function and stability of different proteins, such as the mitochondrial chaperone TRAP1. Interestingly, the SNO site (C501) of TRAP1 is in the proximity of another cysteine (C527). This feature suggests that the S-nitrosylated C501 could engage in a disulfide bridge with C527 in TRAP1, resembling the well-known ability of S-nitrosylated cysteines to resolve in disulfide bridge with vicinal cysteines. We used enhanced sampling simulations and in-vitro biochemical assays to address the structural mechanisms induced by TRAP1 S-nitrosylation. We showed that the SNO site induces conformational changes in the proximal cysteine and favors conformations suitable for disulfide bridge formation. We explored 4172 known S-nitrosylated proteins using high-throughput structural analyses. Furthermore, we used a coarse-grained model for 44 protein targets to account for protein flexibility. This resulted in the identification of up to 1248 proximal cysteines, which could sense the redox state of the SNO site, opening new perspectives on the biological effects of redox switches. In addition, we devised two bioinformatic workflows ( https://github.com/ELELAB/SNO_investigation_pipelines ) to identify proximal or vicinal cysteines for a SNO site with accompanying structural annotations. Finally, we analyzed mutations in tumor suppressors or oncogenes in connection with the conformational switch induced by S-nitrosylation. We classified the variants as neutral, stabilizing, or destabilizing for the propensity to be S-nitrosylated and undergo the population-shift mechanism. The methods applied here provide a comprehensive toolkit for future high-throughput studies of new protein candidates, variant classification, and a rich data source for the research community in the NO field.

Directed Evolution of (R)-2-Hydroxyglutarate Dehydrogenase Improves 2-Oxoadipate Reduction by 2 Orders of Magnitude.

Pathway engineering is commonly employed to improve the production of various metabolites but may incur in bottlenecks due to the low catalytic activity of a particular reaction step. The reduction of 2-oxoadipate to (R)-2-hydroxyadipate is a key reaction in metabolic pathways that exploit 2-oxoadipate conversion via α-reduction to produce adipic acid, an industrially important platform chemical. Here, we engineered (R)-2-hydroxyglutarate dehydrogenase from Acidaminococcus fermentans (Hgdh) with the aim of improving 2-oxoadipate reduction. Using a combination of computational analysis, saturation mutagenesis, and random mutagenesis, three mutant variants with a 100-fold higher catalytic efficiency were obtained. As revealed by rational analysis of the mutations found in the variants, this improvement could be ascribed to a general synergistic effect where mutation A206V played a key role since it boosted the enzyme's activity by 4.8-fold. The Hgdh variants with increased activity toward 2-oxoadipate generated within this study pave the way for the bio-based production of adipic acid.

Unraveling membrane properties at the organelle-level with LipidDyn.

Cellular membranes are formed from different lipids in various amounts and proportions depending on the subcellular localization. The lipid composition of membranes is sensitive to changes in the cellular environment, and its alterations are linked to several diseases. Lipids not only form lipid-lipid interactions but also interact with other biomolecules, including proteins. Molecular dynamics (MD) simulations are a powerful tool to study the properties of cellular membranes and membrane-protein interactions on different timescales and resolutions. Over the last few years, software and hardware for biomolecular simulations have been optimized to routinely run long simulations of large and complex biological systems. On the other hand, high-throughput techniques based on lipidomics provide accurate estimates of the composition of cellular membranes at the level of subcellular compartments. Lipidomic data can be analyzed to design biologically relevant models of membranes for MD simulations. Similar applications easily result in a massive amount of simulation data where the bottleneck becomes the analysis of the data. In this context, we developed LipidDyn, a Python-based pipeline to streamline the analyses of MD simulations of membranes of different compositions. Once the simulations are collected, LipidDyn provides average properties and time series for several membrane properties such as area per lipid, thickness, order parameters, diffusion motions, lipid density, and lipid enrichment/depletion. The calculations exploit parallelization, and the pipeline includes graphical outputs in a publication-ready form. We applied LipidDyn to different case studies to illustrate its potential, including membranes from cellular compartments and transmembrane protein domains. LipidDyn is available free of charge under the GNU General Public License from https://github.com/ELELAB/LipidDyn.

Cancer-related Mutations with Local or Long-range Effects on an Allosteric Loop of p53.

The tumor protein 53 (p53) is involved in transcription-dependent and independent processes. Several p53 variants related to cancer have been found to impact protein stability. Other variants, on the contrary, might have little impact on structural stability and have local or long-range effects on the p53 interactome. Our group previously identified a loop in the DNA binding domain (DBD) of p53 (residues 207-213) which can recruit different interactors. Experimental structures of p53 in complex with other proteins strengthen the importance of this interface for protein-protein interactions. We here characterized with structure-based approaches somatic and germline variants of p53 which could have a marginal effect in terms of stability and act locally or allosterically on the region 207-213 with consequences on the cytosolic functions of this protein. To this goal, we studied 1132 variants in the p53 DBD with structure-based approaches, accounting also for protein dynamics. We focused on variants predicted with marginal effects on structural stability. We then investigated each of these variants for their impact on DNA binding, dimerization of the p53 DBD, and intramolecular contacts with the 207-213 region. Furthermore, we identified variants that could modulate long-range the conformation of the region 207-213 using a coarse-grain model for allostery and all-atom molecular dynamics simulations. Our predictions have been further validated using enhanced sampling methods for 15 variants. The methodologies used in this study could be more broadly applied to other p53 variants or cases where conformational changes of loop regions are essential in the function of disease-related proteins.

DisProt in 2022: improved quality and accessibility of protein intrinsic disorder annotation.

The Database of Intrinsically Disordered Proteins (DisProt, URL: https://disprot.org) is the major repository of manually curated annotations of intrinsically disordered proteins and regions from the literature. We report here recent updates of DisProt version 9, including a restyled web interface, refactored Intrinsically Disordered Proteins Ontology (IDPO), improvements in the curation process and significant content growth of around 30%. Higher quality and consistency of annotations is provided by a newly implemented reviewing process and training of curators. The increased curation capacity is fostered by the integration of DisProt with APICURON, a dedicated resource for the proper attribution and recognition of biocuration efforts. Better interoperability is provided through the adoption of the Minimum Information About Disorder (MIADE) standard, an active collaboration with the Gene Ontology (GO) and Evidence and Conclusion Ontology (ECO) consortia and the support of the ELIXIR infrastructure.

Molecular-dynamics-simulation-guided membrane engineering allows the increase of membrane fatty acid chain length in Saccharomyces cerevisiae.

The use of lignocellulosic-based fermentation media will be a necessary part of the transition to a circular bio-economy. These media contain many inhibitors to microbial growth, including acetic acid. Under industrially relevant conditions, acetic acid enters the cell predominantly through passive diffusion across the plasma membrane. The lipid composition of the membrane determines the rate of uptake of acetic acid, and thicker, more rigid membranes impede passive diffusion. We hypothesized that the elongation of glycerophospholipid fatty acids would lead to thicker and more rigid membranes, reducing the influx of acetic acid. Molecular dynamics simulations were used to predict the changes in membrane properties. Heterologous expression of Arabidopsis thaliana genes fatty acid elongase 1 (FAE1) and glycerol-3-phosphate acyltransferase 5 (GPAT5) increased the average fatty acid chain length. However, this did not lead to a reduction in the net uptake rate of acetic acid. Despite successful strain engineering, the net uptake rate of acetic acid did not decrease. We suggest that changes in the relative abundance of certain membrane lipid headgroups could mitigate the effect of longer fatty acid chains, resulting in a higher net uptake rate of acetic acid.

Ubiquitin Interacting Motifs: Duality Between Structured and Disordered Motifs.

Ubiquitin is a small protein at the heart of many cellular processes, and several different protein domains are known to recognize and bind ubiquitin. A common motif for interaction with ubiquitin is the Ubiquitin Interacting Motif (UIM), characterized by a conserved sequence signature and often found in multi-domain proteins. Multi-domain proteins with intrinsically disordered regions mediate interactions with multiple partners, orchestrating diverse pathways. Short linear motifs for binding are often embedded in these disordered regions and play crucial roles in modulating protein function. In this work, we investigated the structural propensities of UIMs using molecular dynamics simulations and NMR chemical shifts. Despite the structural portrait depicted by X-crystallography of stable helical structures, we show that UIMs feature both helical and intrinsically disordered conformations. Our results shed light on a new class of disordered UIMs. This group is here exemplified by the C-terminal domain of one isoform of ataxin-3 and a group of ubiquitin-specific proteases. Intriguingly, UIMs not only bind ubiquitin. They can be a recruitment point for other interactors, such as parkin and the heat shock protein Hsc70-4. Disordered UIMs can provide versatility and new functions to the client proteins, opening new directions for research on their interactome.

TFG binds LC3C to regulate ULK1 localization and autophagosome formation.

The early secretory pathway and autophagy are two essential and evolutionarily conserved endomembrane processes that are finely interlinked. Although growing evidence suggests that intracellular trafficking is important for autophagosome biogenesis, the molecular regulatory network involved is still not fully defined. In this study, we demonstrate a crucial effect of the COPII vesicle-related protein TFG (Trk-fused gene) on ULK1 puncta number and localization during autophagy induction. This, in turn, affects formation of the isolation membrane, as well as the correct dynamics of association between LC3B and early ATG proteins, leading to the proper formation of both omegasomes and autophagosomes. Consistently, fibroblasts derived from a hereditary spastic paraparesis (HSP) patient carrying mutated TFG (R106C) show defects in both autophagy and ULK1 puncta accumulation. In addition, we demonstrate that TFG activity in autophagy depends on its interaction with the ATG8 protein LC3C through a canonical LIR motif, thereby favouring LC3C-ULK1 binding. Altogether, our results uncover a link between TFG and autophagy and identify TFG as a molecular scaffold linking the early secretion pathway to autophagy.

Characterization of a natural variant of human NDP52 and its functional consequences on mitophagy.

The role of mitophagy, a process that allows the removal of damaged mitochondria from cells, remains unknown in multiple sclerosis (MS), a disease that is found associated with dysfunctional mitochondria. Here we have qualitatively and quantitatively studied the main players in PINK1-mediated mitophagy in peripheral blood mononuclear cells (PBMCs) of patients with relapsing-remitting MS. We found the variant c.491G>A (rs550510, p.G140E) of NDP52, one of the major mitophagy receptor genes, associated with a MS cohort. Through the characterization of this variant, we discovered that the residue 140 of human NDP52 is a crucial modulator of NDP52/LC3C binding, promoting the formation of autophagosomes in order to drive efficient mitophagy. In addition, we found that in the PBMC population, NDP52 is mainly expressed in B cells and by ensuring efficient mitophagy, it is able to limit the production of the proinflammatory cytokine TNF-α following cell stimulation. In sum, our results contribute to a better understanding of the role of NDP52 in mitophagy and underline, for the first time, a possible role of NDP52 in MS.

The conformational and mutational landscape of the ubiquitin-like marker for autophagosome formation in cancer.

Macroautophagy/autophagy is a cellular process to recycle damaged cellular components, and its modulation can be exploited for disease treatments. A key autophagy player is the ubiquitin-like protein MAP1LC3B/LC3B. Mutations and changes in MAP1LC3B expression occur in cancer samples. However, the investigation of the effects of these mutations on MAP1LC3B protein structure is still missing. Despite many LC3B structures that have been solved, a comprehensive study, including dynamics, has not yet been undertaken. To address this knowledge gap, we assessed nine physical models for biomolecular simulations for their capabilities to describe the structural ensemble of MAP1LC3B. With the resulting MAP1LC3B structural ensembles, we characterized the impact of 26 missense mutations from pan-cancer studies with different approaches, and we experimentally validated our prediction for six variants using cellular assays. Our findings shed light on damaging or neutral mutations in MAP1LC3B, providing an atlas of its modifications in cancer. In particular, P32Q mutation was found detrimental for protein stability with a propensity to aggregation. In a broader context, our framework can be applied to assess the pathogenicity of protein mutations or to prioritize variants for experimental studies, allowing to comprehensively account for different aspects that mutational events alter in terms of protein structure and function.Abbreviations: ATG: autophagy-related; Cα: alpha carbon; CG: coarse-grained; CHARMM: Chemistry at Harvard macromolecular mechanics; CONAN: contact analysis; FUNDC1: FUN14 domain containing 1; FYCO1: FYVE and coiled-coil domain containing 1; GABARAP: GABA type A receptor-associated protein; GROMACS: Groningen machine for chemical simulations; HP: hydrophobic pocket; LIR: LC3 interacting region; MAP1LC3B/LC3B microtubule associated protein 1 light chain 3 B; MD: molecular dynamics; OPTN: optineurin; OSF: open software foundation; PE: phosphatidylethanolamine, PLEKHM1: pleckstrin homology domain-containing family M 1; PSN: protein structure network; PTM: post-translational modification; SA: structural alphabet; SLiM: short linear motif; SQSTM1/p62: sequestosome 1; WT: wild-type.

The PyInteraph Workflow for the Study of Interaction Networks From Protein Structural Ensembles.

PyInteraph is a software package designed for the analysis of structural communication from conformational ensembles, such as those derived from in silico simulations, under the formalism of protein structure networks. We demonstrate its usage for the calculation and analysis of intramolecular interaction networks derived from three different types of interactions, as well as with a more general protocol based on distances between centers of mass. We use the xPyder PyMOL plug-in to visualize such networks on the three-dimensional structure of the protein. We showcase our protocol on a molecular dynamics trajectory of the Cyclophilin A wild-type enzyme, a well-studied protein in which different allosteric mechanisms have been investigated.

Structure and Dynamics in the ATG8 Family From Experimental to Computational Techniques.

Autophagy is a conserved and essential intracellular mechanism for the removal of damaged components. Since autophagy deregulation is linked to different kinds of pathologies, it is fundamental to gain knowledge on the fine molecular and structural details related to the core proteins of the autophagy machinery. Among these, the family of human ATG8 proteins plays a central role in recruiting other proteins to the different membrane structures involved in the autophagic pathway. Several experimental structures are available for the members of the ATG8 family alone or in complex with their different biological partners, including disordered regions of proteins containing a short linear motif called LC3 interacting motif. Recently, the first structural details of the interaction of ATG8 proteins with biological membranes came into light. The availability of structural data for human ATG8 proteins has been paving the way for studies on their structure-function-dynamic relationship using biomolecular simulations. Experimental and computational structural biology can help to address several outstanding questions on the mechanism of human ATG8 proteins, including their specificity toward different interactors, their association with membranes, the heterogeneity of their conformational ensemble, and their regulation by post-translational modifications. We here summarize the main results collected so far and discuss the future perspectives within the field and the knowledge gaps. Our review can serve as a roadmap for future structural and dynamics studies of the ATG8 family members in health and disease.

Structure-function investigation of 3-methylaspartate ammonia lyase reveals substrate molecular determinants for the deamination reaction.

The enzymatic reactions leading to the deamination of ß-lysine, lysine, or 2-aminoadipic acid are of great interest for the metabolic conversion of lysine to adipic acid. Enzymes able to carry out these reactions are not known, however ammonia lyases (EC 4.3.1.-) perform deamination on a wide range of substrates. We have studied 3-methylaspartate ammonia lyase (MAL, EC 4.3.1.2) as a potential candidate for protein engineering to enable deamination towards ß-lysine, that we have shown to be a competitive inhibitor of MAL. We have characterized MAL activity, binding and inhibition properties on six different compounds that would allow to define the molecular determinants necessary for MAL to deaminate our substrate of interest. Docking calculations showed that ß-lysine as well as the other compounds investigated could fit spatially into MAL catalytic pocket, although they probably are weak or very transient binders and we identified molecular determinants involved in the binding of the substrate. The hydrophobic interactions formed by the methyl group of 3-methylaspartic acid, together with the presence of the amino group on carbon 2, play an essential role in the appropriate binding of the substrate. The results showed that ß-lysine is able to fit and bind in MAL catalytic pocket and can be potentially converted from inhibitor to substrate of MAL upon enzyme engineering. The characterization of the binding and inhibition properties of the substrates tested here provide the foundation for future and more extensive studies on engineering MAL that could lead to a MAL variant able to catalyse this challenging deamination reaction.

Conformational gating in ammonia lyases.

BACKGROUND: Ammonia lyases are enzymes of industrial and biomedical interest. Knowledge of structure-dynamics-function relationship in ammonia lyases is instrumental for exploiting the potential of these enzymes in industrial or biomedical applications. METHODS: We investigated the conformational changes in the proximity of the catalytic pocket of a 3-methylaspartate ammonia lyase (MAL) as a model system. At this scope, we used microsecond all-atom molecular dynamics simulations, analyzed with dimensionality reduction techniques, as well as in terms of contact networks and correlated motions. RESULTS: We identify two regulatory elements in the MAL structure, i.e., the ß5-α2 loop and the helix-hairpin-loop subdomain. These regulatory elements undergo conformational changes switching from 'occluded' to 'open' states. The rearrangements are coupled to changes in the accessibility of the active site. The ß5-α2 loop and the helix-hairpin-loop subdomain modulate the formation of tunnels from the protein surface to the catalytic site, making the active site more accessible to the substrate when they are in an open state. CONCLUSIONS: Our work pinpoints a sequential mechanism, in which the helix-hairpin-loop subdomain of MAL needs to break a subset of intramolecular interactions first to favor the displacement of the ß5-α2 loop. The coupled conformational changes of these two elements contribute to modulate the accessibility of the catalytic site. GENERAL SIGNIFICANCE: Similar molecular mechanisms can have broad relevance in other ammonia lyases with similar regulatory loops. Our results also imply that it is important to account for protein dynamics in the design of variants of ammonia lyases for industrial and biomedical applications.

S-nitrosylation affects TRAP1 structure and ATPase activity and modulates cell response to apoptotic stimuli.

The mitochondrial chaperone TRAP1 has been involved in several mitochondrial functions, and modulation of its expression/activity has been suggested to play a role in the metabolic reprogramming distinctive of cancer cells. TRAP1 posttranslational modifications, i.e. phosphorylation, can modify its capability to bind to different client proteins and modulate its oncogenic activity. Recently, it has been also demonstrated that TRAP1 is S-nitrosylated at Cys501, a redox modification associated with its degradation via the proteasome. Here we report molecular dynamics simulations of TRAP1, together with analysis of long-range structural communication, providing a model according to which Cys501 S-nitrosylation induces conformational changes to distal sites in the structure of the protein. The modification is also predicted to alter open and closing motions for the chaperone function. By means of colorimetric assays and site directed mutagenesis aimed at generating C501S variant, we also experimentally confirmed that selective S-nitrosylation of Cys501 decreases ATPase activity of recombinant TRAP1. Coherently, C501S mutant was more active and conferred protection to cell death induced by staurosporine. Overall, our results provide the first in silico, in vitro and cellular evidence of the relevance of Cys501 S-nitrosylation in TRAP1 biology.

DisProt: intrinsic protein disorder annotation in 2020.

The Database of Protein Disorder (DisProt, URL: https://disprot.org) provides manually curated annotations of intrinsically disordered proteins from the literature. Here we report recent developments with DisProt (version 8), including the doubling of protein entries, a new disorder ontology, improvements of the annotation format and a completely new website. The website includes a redesigned graphical interface, a better search engine, a clearer API for programmatic access and a new annotation interface that integrates text mining technologies. The new entry format provides a greater flexibility, simplifies maintenance and allows the capture of more information from the literature. The new disorder ontology has been formalized and made interoperable by adopting the OWL format, as well as its structure and term definitions have been improved. The new annotation interface has made the curation process faster and more effective. We recently showed that new DisProt annotations can be effectively used to train and validate disorder predictors. We believe the growth of DisProt will accelerate, contributing to the improvement of function and disorder predictors and therefore to illuminate the 'dark' proteome.

Alterations of the interactome of Bcl-2 proteins in breast cancer at the transcriptional, mutational and structural level.

Apoptosis is an essential defensive mechanism against tumorigenesis. Proteins of the B-cell lymphoma-2 (Bcl-2) family regulate programmed cell death by the mitochondrial apoptosis pathway. In response to intracellular stress, the apoptotic balance is governed by interactions of three distinct subgroups of proteins; the activator/sensitizer BH3 (Bcl-2 homology 3)-only proteins, the pro-survival, and the pro-apoptotic executioner proteins. Changes in expression levels, stability, and functional impairment of pro-survival proteins can lead to an imbalance in tissue homeostasis. Their overexpression or hyperactivation can result in oncogenic effects. Pro-survival Bcl-2 family members carry out their function by binding the BH3 short linear motif of pro-apoptotic proteins in a modular way, creating a complex network of protein-protein interactions. Their dysfunction enables cancer cells to evade cell death. The critical role of Bcl-2 proteins in homeostasis and tumorigenesis, coupled with mounting insight in their structural properties, make them therapeutic targets of interest. A better understanding of gene expression, mutational profile, and molecular mechanisms of pro-survival Bcl-2 proteins in different cancer types, could help to clarify their role in cancer development and may guide advancement in drug discovery. Here, we shed light on the pro-survival Bcl-2 proteins in breast cancer using different bioinformatic approaches, linking -omics with structural data. We analyzed the changes in the expression of the Bcl-2 proteins and their BH3-containing interactors in breast cancer samples. We then studied, at the structural level, a selection of interactions, accounting for effects induced by mutations found in the breast cancer samples. We find two complexes between the up-regulated Bcl2A1 and two down-regulated BH3-only candidates (i.e., Hrk and Nr4a1) as targets associated with reduced apoptosis in breast cancer samples for future experimental validation. Furthermore, we predict L99R, M75R as damaging mutations altering protein stability, and Y120C as a possible allosteric mutation from an exposed surface to the BH3-binding site.

Selective autophagy maintains centrosome integrity and accurate mitosis by turnover of centriolar satellites.

The centrosome is the master orchestrator of mitotic spindle formation and chromosome segregation in animal cells. Centrosome abnormalities are frequently observed in cancer, but little is known of their origin and about pathways affecting centrosome homeostasis. Here we show that autophagy preserves centrosome organization and stability through selective turnover of centriolar satellite components, a process we termed doryphagy. Autophagy targets the satellite organizer PCM1 by interacting with GABARAPs via a C-terminal LIR motif. Accordingly, autophagy deficiency results in accumulation of large abnormal centriolar satellites and a resultant dysregulation of centrosome composition. These alterations have critical impact on centrosome stability and lead to mitotic centrosome fragmentation and unbalanced chromosome segregation. Our findings identify doryphagy as an important centrosome-regulating pathway and bring mechanistic insights to the link between autophagy dysfunction and chromosomal instability. In addition, we highlight the vital role of centriolar satellites in maintaining centrosome integrity.